Efficient algorithms for constraining orientation tensors in Galerkin methods for the Fokker-Planck equation
نویسنده
چکیده
This paper deals with the problem of positivity preservation in numerical algorithms for simulating fiber suspension flows. In contrast to fiber orientation models based on the Advani-Tucker evolution equations for even-order orientation tensors, the probability distribution function of fiber orientation is approximated using the Galerkin discretization of the Fokker-Planck equation with Fourier basis functions or spherical harmonics. This procedure leads to a natural generalization of orientation tensor models replacing ad hoc closure approximations by Galerkin equations for the fine-scale components. After introducing an operator splitting approach to solving the discretized Fokker-Planck equation, we present conditions and correction techniques that guarantee physically correct distribution functions. As the reader will see, the derivation of these conditions is independent of the space dimension and their applicability is not limited to the simulation of fiber suspensions.
منابع مشابه
Pseudo-spectral Matrix and Normalized Grunwald Approximation for Numerical Solution of Time Fractional Fokker-Planck Equation
This paper presents a new numerical method to solve time fractional Fokker-Planck equation. The space dimension is discretized to the Gauss-Lobatto points, then we apply pseudo-spectral successive integration matrix for this dimension. This approach shows that with less number of points, we can approximate the solution with more accuracy. The numerical results of the examples are displayed.
متن کاملAnalysis and Implementation of Numerical Methods for Simulating Dilute Polymeric Fluids
Analysis and Implementation of Numerical Methods for Simulating Dilute Polymeric Fluids David Knezevic Doctor of Philosophy Balliol College Michaelmas Term 2008 In this thesis we develop, analyse and implement a number of numerical methods for simulating dilute polymeric fluids. We use a well-known model in which the polymeric fluid is represented by a suspension of dumbbells in a Newtonian sol...
متن کاملA Fully Discrete Discontinuous Galerkin Method for Nonlinear Fractional Fokker-Planck Equation
The fractional Fokker-Planck equation is often used to characterize anomalous diffusion. In this paper, a fully discrete approximation for the nonlinear spatial fractional Fokker-Planck equation is given, where the discontinuous Galerkin finite element approach is utilized in time domain and the Galerkin finite element approach is utilized in spatial domain. The priori error estimate is derived...
متن کاملNumerical Studies and Simulation of the Lower Hybrid Waves Current Drive by using Fokker – Planck Equation in NSST and HT-7 Tokamaks
Recent experiments on the spherical tokamak have discovered the conditions to create a powerful plasma and ensure easy shaping and amplification of stability, high bootstrap current and confinement energy. The spherical tours (ST) fusion energy development path is complementary to the tokamak burning plasma experiment such as NSTX and higher toroidal beta regimes and improves the design of a po...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Computers & Mathematics with Applications
دوره 71 شماره
صفحات -
تاریخ انتشار 2016